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ABSTRACT 

The processes that form transition disks - disks with depleted inner regions - are not well understood; 
possible scenarios include planet formation, grain growth and photoevaporation. Disks with spatially 
resolved dust holes are rare, but, in general, even less is known about the gas structure. The disk 
surrounding the AO star Oph IRS 48 in the nearby p Ophiuchus region has a 30 AU radius hole 
previously detected in the 18.7 /jm dust continuum and in warm CO in the 5 /jm fundamental ro- 
vibrational band. We present here Submillimeter Array 880 /im continuum imaging resolving an 
inner hole. However, the radius of the hole in the millimeter dust is only 13 AU, significantly smaller 
than measured at other wavelengths. The nesting structure of the disk is counter-intuitive, with 
increasingly large radii rings of emission seen in the millimeter dust (12.9 j^'J AU), 5 pm CO (30 AU) 
and 18.7 /im dust (peaking at 55 AU). We discuss possible explanations for this structure, including 
self-shadowing that cools the disk surface layers, photodissociation of CO, and photoevaporation. 
However, understanding this unusual disk within the stringent multi-wavelength spatial constraints 
will require further observations to search for cold atomic and molecular gas. 

Subject headings: stars: pre-main-sequence — (stars:) planetary systems: protoplanetary disks - 
stars:individual:Oph IRS 48 — planetary systems: formation - submillimeter 



1. INTRODUCTION 

Understanding the interplay between dust and gas is 
crucial for determining the evolutionary pathways of pro- 
toplanetary disks. Gas dominates the mass and dy- 
namics of the disk, but dust is both easier to trace 
observationally and contributes directly to planet core 
formation. Models predict that massive planets clear 
their orbi tal paths, creating distinctive gaps in their na- 
tal disks (iCrida fc Morbidehil [20071: iBrvden et ai1 H99l 
lArtvmowicz fc Lubowl 1 19941 ) . Transition disks, charac- 
terized by observed dust deficits in the i nner regions 
may be the result of this process (e.g. iStrom et al 



1989 1 : iCalvet et al1l2002t iBrown e t al. 2003 iNaiita et al 



20071 ). The recent discoveries of young substellar can- 
didates in the transition disks T Cha and LkCa 15 
lend credence to this hypothesis ([Huelamo et al.l [20111 : 
iKraus fc Irelandl 120121 ) . However, processes other than 
orbital clearing can potentially cause dust deficits, in- 
cluding photoevaporation and grain growth. The gas 
distribution could be diagnostic of different scenarios. 
In disk photoevaporation models, the dust is tied to 
the gas resulting in similar distributions and timescales 
([Alexander et al.1 120061: lOwen et all 120111) . In contrast , 
grain growth to very large sizes would result in an ap- 
parent dust hole but no drop in the gas density. Fi- 
nally, stellar and planetary companions present an in- 
termediate case with some movement of gas across the 
gaps into the inner disk e xpected depending on com- 
panion mass an d location ([Artvmowicz fc Lubowl IT996I : 
IZhu et al.ll20Tl . The presence of gaps can enhance other 
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disk clearin g mechanisms such as the magn eto-rotational 
instability (Chian g" fc Murrav-Clayl I2007D and photoe- 
vaporation ([Alexander fc Armita gc 2009). Knowledge of 
the gas distribution is also particularly relevant for study- 
ing young planets as the gas drives both formation and 
migration. 

The AO star IRS 48 in the p Ophiuchi sta r forma- 
tion region (~ 1 M yr, iLuhman fc Riekd 119991: 121 pc, 
lLoinard et al.l I2008T ) has an unusual disk. Spatially re- 
solved 18.7 pro. imaging of the dust continuum revealed 
a ring-like structure with a peak-to-peak diameter of 
110 AU and a ce ntral hole with a radius of ^30 AU 
(|Geers et al.l [2007). Strong polycyclic aromatic hydro- 
carbon (PAH) emission, commonly a tracer of gaseous 
molecular material, was unresolved and centered on the 
stellar position, interior to the ring seen in the 18.7 fim 
dust. The CO gas might be expected to follow the 
PAH distribution and thus be located inside the dust 
hole. However, resolved imaging with the VLT-CRIRES 
showed warm CO emission in a 30 AU ra dius ring, co- 
incid ent with the proposed dust hole wall ([Brown et al.l 
[20121) . The presence of PAHs was explained by their 
longer survival time relative to CO in high UV flux, low 
density regions. IBrown et al.l ([2012D suggest that trunca- 
tion of the gas and dust, as might be caused by a planet, 
lowers the density in the inner region leading to pho- 
todissociation of CO interior to the truncation radius. 
This scenario should lead to a dust hole detectable with 
millimeter interferometry. 

Millimeter imaging of transitional disks provides some 
of the best constraints on cavity size and characteristics. 
Resolving the cavities requires both long baselines and 
good (u, v) coverage to achieve the necessary high spa- 
tial resolution and dynamic range. The optically thin 
emission observed in the submillimeter is very sensitive 
to the mass surface density profile. However, spatially 
resolved studies have been restricted to the brightest 
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disks. Despite improving statistics, only a limited num- 
ber of transition disks have been observed a t the spatial 
scales necessary to resolve the inner cavity ( Pietu et al 
'2001 iHuehes et al.ll2007l [2001 iBrown et al.l l2008L 1200 
Andrews et alJl200fl l2011bT T~ 



We present in this paper Submillimeter Array (SMA) 
880 /im imaging showing a hole in the distribution of 
large dust grains in the disk around Oph IRS 48. How- 
ever, surprisingly, the size of the millimeter dust hole is 
significantly smaller than that observed in the 18.7 /im 
dust and 5 /im gas. We discuss the observations (jj2J), the 
models used to determine the hole size ($3]) and finally 
discuss the implications of the millimeter data and some 
theories as to the cause of the discrepancy (Q. 

2. OBSERVATIONS AND RESULTS 

We observed Oph IRS 48 (WLY 2-48) with the Submil- 
limeter Array (SMA) in the very extended configuration 
at 345 GHz (880 /im). Two tracks were obtained with 
eight antennas on 2011 August 29 and five antennas on 
2011 August 26. The position measured by the SMA 
was a = 16 h 27 m 37M7, 5 = -24°30'35".5, J2000. Base- 
line lengths ranged from 30 to 590 m. The beam size was 
0'.'45 by 0'.'24 arcsec at a position angle of 27.4°. Double 
sideband (DSB) receivers tuned to 341.165 GHz provided 
4 GHz of bandwidth per sideband. The visibility phases 
were calibrated using the quasars J1625-254 (1.4 Jy) and 
J1626-298 (0.8 Jy), with 10 minute cycles. 3C454.3 was 
used for passband calibration. The absolute flux scale 
was calibrated from a measurement of Neptune from the 
same night. The uncertainty in the flux scale is estimated 
to be 15%. 

Figure [1] presents the image (left) and visibilities 
(right). While the image is clearly not centrally peaked, 
we focus the analysis on the visibilities where a null more 
clearly defines the structure change. Sharp transitions 
in the mass surface density, such as a hole, are seen in 
aperture synthesis observations a s a null in the flux distr i- 
bution with (u,v) distance (e.g.. lAndrews et aLll2011bT ). 
In contrast, most classical T Tauri stars have power- law 
mass surface density profiles, resulting in a smooth de- 
cline of flux with (u,v) distance. Another advantage of 
visibility domain analysis is a better understanding of 
the uncertainties, due to the lack of additional image 
processing, including Fourier transformation and non- 
linear deconvolution. The visibilities in Figure 1 were 
deprojected to account for position angle, PA, and disk 
inclination, i, and b inned assuming a circular disk (e.g., 
iHughes et alJ 12007). The deprojected (u,v) distance is 
R = V^a + w ith d& = Rsin<f> and db = Rcos(f>cosi 
where <fi — arctan(v/u) - PA. 

Earlier archival data from 2007 May 14 observed IRS 
48 at 230 GHz (1.3 mm) in the very extended config- 
uration. The track had only 6 working antennas and a 
beam size of 0767x074 (P.A. = -7.8°). The spatial resolu- 
tion was inadequate to clearly resolve the millimeter dust 
hole, but the data indicated that the millimeter dust was 
unlikely to follow the shorter wavelength emission. The 
visibilities are overlaid in gray on Figure [1] but all further 
analysis uses the higher spatial resolution 880 /im data. 

Spatial filtering of flux from the outer disk beyond 
~100 AU is a concern with the long baselines in the SMA 
data. IRS 48 has been observed with single dish millime- 



ter telescopes with a flux of 180±9 mJy with JC MT- 
SCUBA at 860 /im (lAndrews fc Williams! f2007h and 
60±1 mJy with the IRAM 30m at 1.3 mm (jMotte et all 
Il998h . The SMA tracks have total fluxes of 160 mJy at 
880 /xm and 50 mJy at 1.3 mm. Less than 20% of the 
total flux has thus been lost in the interferometric data, 
ruling out the possibility of a large reservoir of outer disk 
millimeter dust. 

3. MODELING 

We model the millimeter emission following the simi- 
larity solution approach which has been used to charac- 
terize the millimeter dust in both normal and transit ion 
disks (IHughes et al.|[2008t lAndrews et al.ll20"09ll2011allbh . 
Our emission model is taken to be azimuthally symmetric 
with a surface brightness profile appropriate for optically 
thin thermal emission, 



I u k B„{T d )(l - e~ T ) ^B v (T d ) T 



(1) 



where B u is the Planck function, T d is the dust temper- 
ature, and r is the optical depth. The dust temperature 
is assumed to follow a power-law profile normalized by 
the stellar temperature of IRS 48: 



-0.5 



(2) 



where T»=9000 K and r*=1.35 R . The optical 
depth across the disk is set by the surface density for 
which we adopt the simil arity solution appropriate for 
a viscous accretion d isk (|Lvnden-Bell fc Pringld Il974t 
lHartmann et al.lll998[ ): 



r cx 



exp 



2~7 



(3) 



where 7 sets the gradient of the profile and r c is the char- 
acteristic scaling radius. We include a cavity by setting 
to zero any emission within its radius, i? C av 

After normalizing and projecting the model image onto 
the sky plane with some inclination, position angle, and 
positional offsets, we use MIRIAD to sample the model 
visibilities. The position angle is fixed to 90° based 
on the shorter wavelength data. The fitness of each 
model, or the likelihood, is constructed from the chi- 
squared statistic assuming the visibilities have uncorre- 
cted, Gaussian errors, L oc exp(— x 2 /2). To estimate 
the best fit model parameters and their associated un- 
certainties, we sample the pos terior distribution us i ng an 
ensemble MCMC algorithm ()Goodman fc Weard l2010l: 
Forcman-Mack ev et al.ll2012D . We report the best fit pa- 
rameters and the uncertainties from the 68% confidence 
intervals of the marginalized distribution (Table [T]). The 
best fit model image and visibility fits are shown in Fig. 

IB 

We also explored a couple of alternate models to un- 
derstand the constraints the data put on the millimeter 
dust distribution. First, a constant density ring model 
required flux between 7 and 48 AU to match the visibili- 
ties. Second, R rav was fixed at 3 AU to match the warm 
CO gas ring (|Brown et al. 2012) and a constant density 
of dust was introduced interior to i? C av This model still 
required 50% of the millimeter flux to come from within 
i? cav . In both cases, the fit produced looked adequate by 
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Fig. 1. — Images of the data (left), model (center left) and residuals (center right). The contours are 3<r intervals where the RMS noise 
is 2.7 mjy/beam. The (u,v) visibities for 345 GHz (black) and 230 GHz (gray) are on the far right with the best fit model from Tablc[l]m 
red. 



TABLE 1 

Model parameters 



Parameter 


Best fist 


Lower bound 


Upper bound 


Fsso^m (Jy) 


0.16 


0.150 


0.165 


r c (AU) 


35 


27 


36 


7 


0.12 


-0.49 


0.6 


H cav (AU) 


12.9 


9.5 


14.6 


Inclination (°) 


35 


28 


38 



eye but had a larger \ 2 than the similarity solution fit 
reported in Table [TJ The implication is that, regardless 
of model details, both a central dropoff in the millimeter 
emission and substantial dust within the 5 /an CO ring 
are required. 

Both the 18.7 /im and CO images show asymmetries 
with the north side brighter than the south. However, 
the millimeter image shows no significant asymmetries 
beyond those expected from an inclined disk. The asym- 
metries at shorter wavelengths may reflect variable ex- 
tinction of the illuminating radiation field. 

The derived inclination of 35° is smaller than those 
calculated from the 5 fim CO (42°±6°) and the 18.7 /im 
dust (48° ±8°). The inclination was left as a free param- 
eter in the similarity solution models with lower inclina- 
tions resulting in a smoother visibility profile and smaller 
X 2 ■ This could be the result of azimuthal averaging over 
asymmetric structure that is unseen at the SMA resolu- 
tion. Large inclinations result in a slightly larger central 
cavity in the millimeter dust with 45° corresponding to 
a best fit cavity size of 15 AU. 

4. DISCUSSION 

Oph IRS 48 is well-studied with resolved measurements 
of the disk and central hole at several wavelengths and, 
unusually, constraints on both the gas and dust. Despite 
spatially resol yed rings of emiss ion in the 18.7 /im dust 
(55 AU p eak: IGeers et alJ 12001 and the 5 /im CO (30 
AU peak; iBrown et alj|2012t ). the inner edge of the mil- 
limeter dust lies interior to both, with a cutoff in the dust 
density at the even smaller radius of 12.9 ±\'l AU (Fig. 
[2|). The SMA image definitively rules out the millimeter 
emission following the structure in either the 5 /im CO 
or the 18.7 /im dust. The central hole is barely resolved 
with the SMA with a large portion of flux within the 
central beam. Thus, all models require a significant per- 
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Fig. 2. — The radial distribution of flux from the millimeter dust 
(black), 5 fim CO (red) and 18.7 fim dust (blue). The millimeter 
profile is the best fit profile listed in Table [l] with the grey region 
representing fits wit hin the error bars. T he 5 fim CO profile is the 
best fit model from IBrown et all 1120121) Figu re 8. The 18.7 fim 
dust flux is a cut along the major axis of the [Gccrs et al.] 1120070 
VISIR image after 2-D maximum likelihood deconvolution of the 
PSF. 

centage of the millimeter flux to come from within the 
CO and IR dust rings. If the disk is currently forming 
planets, the complex disk structure may provide regions 
advantageous to planet and planetesimal formation. It 
may be the case that already formed planets are respon- 
sible for some of the sculpting seen in the disk. 

Both infrared measurements likely trace only the sur- 
face layers of the disk. The VLT-CRIRES 5 /im CO 
observations are very sensitive to warm CO with con- 
straints both spatially and spectrally However, the data 
are less sensitive to cold CO and completely insensitive to 
non-molecular gas or other species of molecular gas. The 
VISIR 18.7 /im image has a PSF of ~0 / /5 which is similar 
to the CO hole size. The 18.7 /tm dust hole may be larger 
than the 30 AU previously reported ([Gecrs "etall 120071) 
as seen in the deconvolved cut in Figure [5] (noted also 
by Maaskant et al., in prep; Bruderer et al., 2011, priv. 
comm.). While the majority of the 18.7 /im emission is 
from outside the millimeter dust, the inner regions may 
not be entirely void of infrared emitting dust. 

The millimeter dust emission places some constraints 
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on the amount of gas within the hole under assumptions 
of a canonical gas to dust ratio. In the best fit model, 
71% of the millimeter flux arises from dust within the 30 
AU radius gas ring. From all the models, 40-90 % of the 
mil limeter emission must be interior to 30 AU. We use 
the iBeckwith etall $1990) conversion of millimeter flux 
to disk mass with a f3 of 1 (jRodmann et al J 12006). The 
dust temperature is taken to be in the range of 50-150 K. 
The gas mass within the 30 AU CO ring is then 10 -4 - 
10 -3 M Q corresponding to a CO mass of 10~ 7 -10~ 6 M , 
assuming a H2/CO number ratio of 10 4 . In a flared disk 
with a scale height set to the hydrostatic pressure at 200 
K, the density of CO is 10 6 -10 7 cm -3 or a surface density 
of 10 19 -10 20 cm -2 , comparable to a disk with no cavity. 
Significantly warmer gas is likely to produce 5 /im CO 
emission detectable in the CRIRES data. 

There are three basic scenarios to reconcile the mil- 
limeter dust and 5 /zm CO emission: 1.) the gas (and 
dust) is cold and therefore not detectable in the infrared, 
2.) the CO is photodissociated and only atomic gas is 
left interior to 30 AU, or 3.) the gas has been removed 
leaving millimeter dust and PAHs. 

4.1. Scenario 1 

Shadowing from a puffed up inner rim, either at the 
millimeter dust wall or from atomic gas or PAHs closer to 
the star, could cool the surface layers of the disk further 
out. The 5 /im CO limits the amount of warm gas within 
the 30 AU gas radius but the l imits are much less strin- 
gent at cooler temperatures (jBrown et al.1 120121 their 
Figure 12). A temperature of ^150 K or less through- 
out the inner disk is needed to reconcile the 5 /jm CO 
surface density limit with the modeled millimeter sur- 
face density, assuming a gas-to-dust ratio of 100. Any 
illuminated wall would also likely need to be devoid of 
CO. The shadowing region would require a large scale 
height to shadow the disk out to 30 AU. In order for this 
scenario to be plausible, the inner regions must be very 
settled while the outer structure is flared, exposing the 
surface layers to stellar radiation to produce the infrared 
emission. 

4.2. Scenario 2 

Strong enough UV fields can photodissociate CO, with 
the photodissociation rate in an unattenuated interstel- 
lar field being 2 x 10" 10 s -1 . The UV field around an AO 
star such as IRS 48 is much stronger than the interstel- 
lar field even in the critical 912-1100 A range and would 
lead to a photodissociation lifetime of <1 year. However, 
CO photodissociation is subject to self-shielding with 
the lines becoming saturated at a 12 CO column density 
of about 10 15 cm" 2 , strongly decreasing the photodis- 
socia tion rate (|van Dishoeck fc "B lack 1988; Vissc r et al.l 
2009). Both millimeter dust grains and PAH molecules 
are also effective absorbers of UV radiation, lowering the 
UV field strength and thus the photodissociation rate. 
The estimated CO density based on the millimeter dust 
would be strongly self-shielding making complete pho- 
todissociation unlikely. 

4.3. Scenario 3 

Removal of only the gas within 30 AU resulting in a 
strongly decreased gas-to-dust ratio is a final possibility. 



Photoevaporation primarily affects the gas with small 
dust grains tied to the gas flow. It may be possible 
for large enough dust grains to remain. Destructive 
collisions of the remaining large particles would then be 
needed to generate the PAHs also seen from the inner 
regions. 



These different scenarios should leave differing obser- 
vational traces which further observations may be able to 
distinguish. In the case of cold CO, ALMA observations 
will have the sensitivity and spatial resolution to look 
for cold CO within the central region of the disk. For 
photodissociated CO, optical spectroastrometric or high 
spectral resolution submillimeter observations of atomic 
C or O would determine if the photodissociation products 
of CO are present within 30 AU. If the gas is removed 
completely, neither cold CO nor the atomic components 
of CO would be present in the inner disk region. 

The ring of 18.7 /an dust is even more difficult to un- 
derstand with most of the emission coming from outside 
either the millimeter dust or the 5 /im CO. The dominant 
18.7 /im flux from ~55 AU has no obvious counterpart 
in the millimeter emission. Like the 5 /im CO, the 18.7 
/im flux comes from the disk surface and could reflect dif- 
ferences in illumination and heating of the surface layers 
rather than underlying changes in the distribution. How- 
ever, the centrally-peaked PAH emission points to some 
heating and small particles in the inner region. Recent 
studies in the mid-IR and scattered light suggest that 
small surface partic les within the cavities seen in large r 
grains are common (|Muto et alll2(]IllDoriFeraIl l2012). 
Holes visible in the mid-IR are often of similar size to 
those seen in the sub millimeter (e.g. iPanic et al.l 120121 : 
iThalmann et al.ll2010l ). However, the scattered light pro- 
files, tracing the small surface grains, are smooth with 
no decrease within the cavity. Dust filtration or trap- 
ping creating grain size differences across planetary gap s 
has been proposed (|Rice et al.ll2006t iPinilla et al.lf2012fl . 
with the largest particles left on the outside of the gap 
and smaller particles, being more closely tied to the gas, 
allowed closer to the star. While the mid-IR dust and 
PAHs in Oph IRS 48 may fit this paradigm, the large 
millimeter dust interior to the smaller mid-IR emitting 
dust does not. Likewise, 5 /im CO emission usually arises 
interior to s ubmillimeter holes rather than from furthe r 
out in disk (jSalvk et al.l 120091 : IPoiitoppidan et al.ll2008f ). 
Continued multiwavelength studies are vital to search for 
other examples of the inside-out nesting structure seen 
in the Oph IRS 48 disk. 

5. CONCLUSIONS 

We have used new spatially resolved SMA 880 /im ob- 
servations to place stringent constraints of the dust dis- 
tribution within the Oph IRS 48 disk. The SMA data 
clearly rules out that the millimeter dust distribution fol- 
lowing either the 5 /im CO or 18.7 /im dust . A large 
fraction of the millimeter flux is from dust interior to 
both the 18.7 /im dust ring (55 AU radius) and the 5 
/im CO emission (30 AU radius). The millimeter dust 
also shows a strong drop in surface density with a null 
in the visibilities. Our similarity solution models place 
the hole at 12.9 I34 AU. The millimeter emission clearly 
shows that an empty cavity is not responsible for the CO 
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gas ring described in lBrown et al.l (|2012l ). Further obser- 
vations are needed to distinguish between the potential 
scenarios of cold shadowed gas, photodissociated CO and 
low gas-to-dust ratio, but all these scenarios result in an 
unusual environment for potential planet formation. 

The Submillimeter Array is a joint project between 
the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astro- 
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the Smithsonian Astrophysical Observatory for support 
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